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Abstract 

<D The predictions of Drell-Yan production of low-mass, lepton-pairs, at high rapidity 

at the LHC, are known to depend sensitively on the choice of factorization and renor- 
^ malization scales. We show how this sensitivity can be greatly reduced by fixing the 

^ factorization scale of the LO contribution based on the known NLO matrix element, so 

that observations of this process at the LHC can make direct measurements of parton 
distribution functions in the low x domain; x 5~ 10~^. 
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> 1 Introduction 

X 

^ The very high energy of the LHC allows us to probe the parton distribution functions (PDFs) 
of the proton at extremely small x, a region not accessible at previous accelerators. To extract 
the PDFs we describe the experimentally observed cross sections as a convolution of the PDFs 
and the cross section for the hard subprocess. Schematically, we may write the cross section in 
the form 

da/d^p = J dxidx2 PDF{xi,fiF) \Mip;iiF,I^R)\^ PDF(x2,/xf) , (1) 

where a sum over the various pairs of PDFs is implied. The matrix elements squared, 
describe the cross sections of the elementary partonic subprocesses. An observable process at 
the LHC samples PDFs that carry momenta fractions Xi of that of the initial protons, where 

a;i,2 = ^exp(±r). (2) 
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Here mhard denotes the mass created in the subprocess. We see that if a relatively low-mass 
subsystem is observed with large rapidity, Y, at large collider energy, then we probe PDFs at 
very small values of X2- However, the problem is that using ([T| we do not know the factorization 
scale fip which the PDFs are measured. Moreover, in the low x region, the PDFs strongly 
depend on the choice of /xj^Q It is made worse due to the dominance of the gluon PDF at 
small X. Consider, for example, the Drell-Yan production of a low-mass, M, lepton pair. For 
a relatively low factorization scale the LO qq — )■ 7* subprocess may be overshadowed by the 
NLO subprocess gq — )■ 57*, due to the dominance of the gluon PDF at low x. However, as we 
shall see below, it is this very dominance which will allow us to introduce a procedure which 
greatly suppresses the scale dependence of the predictions. 

In general, after the summation of all perturbative orders, the final result should not depend 
on the choice of fip that is used to separate the incoming PDFs from the hard matrix element. 
Contributions with low virtuality, < of the incoming partons are included in the PDFs, 
while those with > fi'jp are assigned to the matrix element. However, at low x the probability 
to emit a new parton in some Ahf interval is enhanced by the large value of the longitudinal 
phase space, that is by the large value of ln(l/a;). In fact, the mean number of partons in the 
interval Alufip is 

{n) ~ ln(l/x) Aln/i|, (3) 

71 

leading to a value of (n) up to about 8, for the case ln(l/x) ~ 8 and the usual fip scale variation 
from yu/2 to 2/i. On the other hand, the NLO coefficient function (the hard matrix element) 
allows the emission of only one parton. Therefore we cannot expect compensation between the 
contributions coming from the PDF and the coefficient function. (At large x the compensation 
is much more complete and provides reasonable stability of the predictions to variations of the 
scale fip-) 

Here, our plan is to use the NLO result to fix the proper choice of the factorization scale 
for the LO part of the amplitude, and to demonstrate that the resulting predictions are quite 
stable to variations of fip in the remaining NLO part. In other words, we would like to choose 
the value of fip which minimizes the higher order NLO, NNLO,.. contributions. 

As mentioned above, we may probe PDFs at low x by observing low-mass systems, such 
as Drell-Yan lepton-pair production or bb production, at the LHC. In each case the subprocess 
mass in ^ is relatively small, and so observation at large rapidities will probe PDFs at small 
values of X2. These kinematics are particularly optimal for the LHCb experiment |3]. Since 
Drell-Yan is the simpler example, we will use it to illustrate the method. 

To sketch the idea, we start with the LO expression for the cross section. In the coUinear 
approach, the cross section has the form 

aifip) = PDF(/i^)®C^^®PDF(/ii.). (4) 

^See the dependence of low-mass Drell-Yan production or bb production, shown, for example, in [P and 
|2] respectively. 
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The effect of varying the scale from mto/j^F, in both the left and right PDFs, can be expressed, 
to first order in a^, as 



aifxp) = PDF(m) ® ( C 
where the splitting functions 
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act on the right and left PDFs respectively, see Fig. [TJ We may equally well have incoming 
g's in Pright and incoming g's in Pieft. Recall that in calculating the correction in (|5]), the 
integral over the transverse momentum (virtuality) of the parton in the LO DGLAP evolution 
was approximated by the pure logarithmic dk'^/k'^ form. That is, in the coUinear approach, the 
Leading Log Approximation (LLA) is used. 
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(a) qq ^ gy* (b) gq ^ qy* (c) NNLO 

Figure 1: Diagrams (a,b) show NLO subprocesses for Drell-Yan production resulting from the 
splitting of the upper PDF (or 'right' PDF, in the notation of ([5])). Diagram (c) is the main NNLO 
subprocess gg — )■ g7*g. 



Now let us study the expression for the cross section at NLO. First, we note that the 
original Feynman integrals corresponding to the NLO matrix element, C^^^, do not depend on 
fj,F- However, we will see below how a scale dependence enters. At NLO we may write 

a{^IF) = PDF(^^)® (C^o + a,C,^„^0)®PDF(/ii.), (7) 

where we include the NLO correction to the coefficient function. In terms, for example, of the 
diagrams (a,b) of Fig. [T] this means that the 2 — )■ 2 subprocesses qq — )■ gj* and gq — )■ q'j* are 
now calculated with better, than LLA, accuracy. However part of this contribution was already 
included, to LLA accuracy, in the second term in ([s]). So this part should be subtracted 
from C^^'^. Moreover, this LLA part depends on the scale fip- As a result, changing fip 
redistributes the order correction between the LO part (PDF (g) C^*"* CS> PDF) and the NLO 
part (PDF ® agC^^]^ (g) PDF). We see the part of NLO correction which remains after the 
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subtraction, C^f^ifip), now depends on the scale fip, coming from the fip dependence of the 
LO LLA term that has been subtracted off. The trick is to choose an appropriate scale, fip = fJ'O, 
so to minimize the remaining; NLO contribution C^^^f^ifip)- To be more precise, we choose a 
value hf = 1^0 such that as much as possible of the 'real' NLO contribution (which has a 
'ladder-like' form and which is strongly enhanced by the large value of ln(l/x)) is included in 
the LO part where all the logarithmically enhanced, asln(l/a:), terms are naturally collected 
by the incoming parton distribution^ 

In the next section we describe how to choose an appropriate scale Hf = A^o- Then in 
Section |3] we demonstrate how by fixing /li? = /xq in the LO part of (|5|, we obtain stabihty of 
the Drell-Yan cross section with respect to variations of ^f the remaining NLO contribution. 
Moreover, it turns out that the NNLO correction to the Drell-Yan process at such a value of 
/ii? is also very small. 



2 The choice of an appropriate factorization scale 

DGLAP evolution describes the variation of parton densities a{x,fi) with the scale /x. At LO 
it accounts for the splitting of a parton a{xa, kat) to a parton h{xh, ht), with the longitudinal 
momentum fractions which satisfy, Xa > Xb and the transverse momenta which satisfy kat < kbt- 
This splitting is described by Feynman diagrams in the Leading Log Approximation (LLA) 
where the exact kt distribution is approximated by the pure logarithmic form dkf/k^. 

However the hard matrix element is calculated more precisely. Here we account for the exact 
kt dependence. Our aim is to choose the upper limit /i^? in the logarithmic kf integration in 
such a way that the value of the integral becomes equal to the result given by the NLO matrix 
element which contains the exact kt integration. Strictly speaking this cannot be done for all 
the hard subprocesses simultaneously, with a common value of fiF- So we have to compromise, 
and to choose ^f for the subprocess which gives the major contribution. 

For Drell-Yan production at low x, the majority of quarks/antiquarks are produced via the 
low-x gluon to quark splitting, Pgg. That is, the most important NLO subprocess is gq — )■ 57*. 
The corresponding cross section reads 

where = 3 and u = — s — t. The kinematical upper limit of \t\ is 

|t| = s - = (1 - z)MV^, (9) 



^Actually our approach is rather close to the fct-factorization method. Using the known NLO result we 
account for the exact kt integration in the last cell, adjacent the LO 'hard' matrix element (describing heavy 
photon emission), while the unintegrated parton distribution is generated by the last step of the DGLAP 
evolution, just like the prescription proposed in [4ji5j. 
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where M is the mass of Drell-Yan lepton pair ancj^z = M'^/s. 

At first sight, a problem appears to be that, for a large incoming gq energy, V^, we will 
have to consider very large \t\; that is unreasonably large values of /e^ = — (1 — z)t. However, 
in reality this is not the case. Let us consider the previous splitting, Pgg, which produces the 
gluon from an incoming quark. That is the qq — > qqj* subprocess. For large energy, Sqq ^ M^, 
of the incoming qq system, the gluon flux generated by the fast quark is of the form dxg/xg. 
Since the cross section ([s]) decreases as the subenergy s increases, weighting the gq — )■ 57* cross 
section with the gluon flux, results in the main contribution coming from the relatively low 
s ~ region. 

An explicit calculation shows that in such a case the pure logarithmic LO evolution inte- 
grated up to yUi? = /xo = 1.4M convoluted with C^'^ reproduces the NLO contribution of this 
subprocess. Since in the low-x region gq — )■ 57* is the dominant subprocess we anticipate that 
the choice of /i^? = /iq = 1.4M will minimize the remaining NLO contribution. 

Moreover, at NNLO, for very small x, the dominant diagram is gg — > q'~f*q, where on both 
sides of 7* (the Drell-Yan heavy photon) the quark and the antiquark are produced via the 
gluon to quark splitting, see Fig. [l[c). That is, we have a gg ^ 11*1 subprocess with the same 
kinematics as for the NLO qg — ?■ g7* case considered above. It should be also well mimicked 
by the LO evolution with fip = /iQ. Thus the choice fip = P-o will suppress (in a low-x domain) 
the NNLO contribution as well. 



3 Scale dependence of low-mass Drell-Yan production 

Figs. [2] and |3] show, respectively, the factorization, fip, and the renormalization, hji, scale 
dependences of the cross section as a function of rapidity for the Drell-Yan productioi]|^ of a 
fi~^fi~ pair of mass M = 6 GeV at the LHC energy of 7 TeV. As the plots are symmetrical 
about rapidity F = 0, we are able to compare directly two sets of distributions. 

The left-half of each of these plots shows the LO, NLO and NNLO results calculated with the 
MSTW2008 LO, NLO and the NNLO PDFs [7] respectively. In Fig. ga) the renormahzation 
scale = M is fixed, while the value of factorization scale is varied: fip = M/2, M, 2M. 
In Fig. |3](a) the factorization scale is fixed fip = M, while the renormalization scale is varied: 
Hr = M/2, M, 2M. It is clear that the uncertainties arising from the fip dependence are huge 
at the LO, large at the NLO and still not negligible even at the NNLO; see, also, [1]. The 
dependence on the value of fiR is not so strong, but is still not negligible. 

In the right-half of these plots, the continuous bold lines show the cross sections calculated 
using the PDFs with our 'optimal' choice of scale, fip = fiQ = 1AM, in the LO part. The only 

^Strictly speaking, z is the ratio of the hght-cone momentum fraction carried by the 'daughter' quark to that 
carried by the 'parent' gluon, z = x+/a;+. 

■*We use the computer code Vrap of Ref. [6] . 
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Figure 2: (a) Sensitivity of M = 6 GeV Drell-Yan production, as a function of rapidity Y, 

to the choice of factorization scale: fip = M/2, M, 2M, at LO, NLO, NNLO. (b) The bold lines 
correspond to the choice jj,p = Hq = 1AM which minimizes Cf^, and show the stability with 
respect to the variations /i = M/2, M, 2M in the scale of the PDFs convoluted with C^^^o _ 
11 dependence is indicated by the symbolic equation at the top of the diagram. The dashed lines 
show that the stability disappears for other choices of jiQ. The small crosses are the NNLO result. 
In this figure the renormalization scale is fixed at ijlr = M. 



scale dependence in Fig. comes from the variation of fip = M/2, M, 2M in the PDFs 
convoluted with the remaining NLO contribution. Recall that the value of /io was found to 
almost nullify the NLO correction, leaving only a small C^^]f contribution. Indeed, the plot 
shows that the distribution has very good stability to the variation of fip. The dashed lines on 
this plot (Fig. ^h)) show that the stability would be much worse if we were to take /xq = M 
or fiQ = 2M, rather than the optimal value of fip = fJ'O = 1.4M. 

Moreover, as we argued in Section |2| even the NNLO correction is rather small for the choice 
IIq = 1AM. The corresponding NNLO distribution, calculated with fip = fiQ = 1AM, is shown 
by the crosses. It is rather close to the (NLO) continuous lines. Recall, here we have used 

= M, and the same MSTW2008 NLO PDFs for both the NLO and NNLO calculations, as 
we wish to see the effect of changing from C^^'~' to C*'^^'"^ leaving the PDFs unchanged. Let us 
explain why this is the appropriate procedure, since at first sight it may appear that we should 
use NNLO PDFs for the NNLO prediction. 
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Figure 3: (a) Sensitivity of M = 6 GeV Drell-Yan /x+yU production to the choice of renormalization 
scale: hr = M/2, M, 2M, at LO, NLO, NNLO; and the fixed factorization scale /xp = M. (b) 
Stability of the NLO result is achieved for the optimal choice fip = l^o = 1.4M. 



The basis of our method is that the long evolution length, in ([5]) with m = Qo, means that 
the LLA LO result can capture most of the answer provided that the factorization scale is chosen 
correctly. That is, the LO term captures most of the NLO contribution ior fip = fiQ = 1.4M, 
and even the major part of the NNLO contribution for the same value of fip. 

In general, one may use the known NNLO result to determine another value of the optimal 
scale. Hi, such that the NLO part (the term in the splitting functions) of DGLAP evolution 
will capture most of the NNLO contribution. However, as it is seen from Figs. [2]and|3| already 
the scale /iq used for the LO part (the a° term in the splitting functions) provides sufficient 
stability of the results for the Drell-Yan cross sectioiiQ 

In Fig. |3](b) we see an analogous stability with respect to variations of the renormalization 
scale: /ir = M/2, M, 2M, when the factorization scale is fixed at fip = f^o = 1.4M. This 
additionally confirms that the NLO (and NNLO) corrections are suppressed for the optimal 
choice of fip. 

^Recall, that the subtractions of the contribution generated by the NLO DGLAP evolution convoluted with 
the LO matrix element from the exact NNLO matrix element include the NLO splitting functions. Therefore 
we can convolute the known NNLO coefficient functions, C'NNLO^ with either the NLO or NNLO PDFs, but not 
with the LO PDF; in the last case we need to recalculate the NNLO matrix element. 
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4 PDF dependence of low-mass Drell-Yan production 



What values of Xi and X2 of the PDFs are probed in Drell-Yan production of a system of mass 
M with rapidity Y7 Formally, in the coUinear approximation, the LO term samples PDFs with 
precisely known momentum fractions 



M 



exp(±r) . 



(10) 



That is, the Xi distributions are 5-functions. This is not true in reality. Formula (10) does not 



include the transverse momenta of the incoming partons. It assumes kt -C M. If we account 
for kt, then we would have 

v/M2 + k1 



Xl,2 



exp(±r) . 



(11) 



On the other hand, the optimal scale no = 1AM is rather large. For this scale, the partons in 
the PDFs have kts up to /xq, which may exceed the value of M. In order to get a more realistic 
impression of the interval of x sampled by the process, we go back one step in the DGLAP 
evolution. That is, we consider the qg ^ subprocess, which we already discussed at NLO. 
Recall that actually the majority of low-x quarks(antiquarks) were produced just in this way. 
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Figure 4: The values of xi and X2 of PDFs probed by the observation of Drell-Yan production, 
through the qg channel, of a lepton-pair of mass M = 6 GeV at a 7 TeV LHC. The gluon (quark) 
PDF is sampled in the x region given by the dashed (continuous) curves. The distributions are 
shown for two different values of the rapidity, Y , of the lepton pair system. The two curves for each 
PDF at each rapidity correspond to two choices of the virtuality cutoff: namely Qq = 1 and Qq = 2 
GeV. 
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To be specific, let us take the gq — )■ 57* matrix element (|8j), and convolute it with the NLO 
MSTW PDFs. Fig. |4] shows the resulting distributions of the momentum fractions of the quark 
(continuous lines) and the gluon (dashed lines). As before we consider Drell-Yan production of 
a system of mass M = 6 GeV at the LHC energy ^/s = 7 TeV. We show the distributions for 
two rapidities: Y = and Y = 4. The very sharp left-hand cutoff of the distributions refiects 



the fact that the momentum fraction cannot be less than that given by (10). The tail at larger 
X is due to the smearing of the distribution arising from the {qg — )■ 7*g) subprocess. A narrow 
peak near the lowest possible value of x is due to the 1/t singularity in the subprocess cross 
section 



To calculate the distributions shown in Fig. |4]we introduced the infrared cutoff |t| > Qq, 
with Qq = 1 GeV (the smallest virtuality at which the DGLAP evolution starts). For a larger 
cutoflQ Qo = 2 GeV we get the lower peaks shown in the diagram. For the cutoff Qq = 1 GeV, 
the subprocess gq — )■ 57* gives more than 90% of the whole Drell-Yan cross section. Thus the 
distributions shown in Fig. |4] provide a rather good evaluation of the x-regions sampled by the 
Drell-Yan process for these kinematics. 

What information about PDFs do we obtain from measurements of low-mass Drell-Yan 
lepton-pairs at high rapidity? It is important to emphasize that the data make a direct mea- 
surement of the quark/ ant iquark PDFs at the optimal scale fi = 1.4M, with little factorization 



or renormalization scale ambiguity, for the values of Xi^2 given by (10). That is, there is a direct 
one-to-one correspondence between the data and the PDFs. For example, the production of 
a low-mass M = 6 GeV pair at 7 TeV, going forward with rapidity Y = 4, measures 

quark/antiquark PDFs with 

xi = 4.7 X 10"^ X2 = 1.6 X 10-^ (12) 

but at a moderately high scale, /i^ ~ 70 GeV^. 

What is the impact of such measurements on global PDF analyses? Clearly, the measure- 
ments can probe a small x domain so far unexplored by data. If, at such small x, we were to 
believe in pure DGLAP evolution, then the majority of the low x quarks (antiquarks) comes 
from the g ^ qq splitting. Therefore the M dependence of the Drell-Yan cross section will 
probe the gluon distribution in a similar way to how the scaling violations in deep inelastic 
scattering data, dF2{x, Q^)/dQ'^, probe, via DGLAP evolution, the low x gluon densitjj^ So, by 
including low-mass high-rapidity Drell-Yan data in the global parton analysis we will strongly 
constrain the low x gluon PDF. Note that Drell-Yan data at the LHC may probe very low 
values of x. In this future global analysis we should account for the absorptive, ln(l/x), etc., 
modifications to the evolution equations, as well as using the Drell-Yan data at the optimal 
scale. 



^In the case of Qo we come back to (5-functions at the x values given by (10). 
^Recall, that starting the evolution from Qq — 1 GeV more than 90% of the 
described by the gq — j*q subprocess, and so will probe g{x2)q{xi), see Fig. 4. 



Drell-Yan cross section is 
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Figure 5: The 'LO+NLO' cross sections for the Drell-Yan production of a Af = 6 and M = 12 GeV 
pair, as a function of its rapidity, at 7 TeV, obtained using four different recent NLO sets of 
PDFs in El El [ID]- We show the la error corridor for the predictions obtained using the MSTW 
parton set. 



Bearing in mind the above qualifications, let us use recent sets of existing global PDFs to 
calculate the cross section for the examples of M = 6 and M = 12 GeV lepton-pair production. 
The results are shown in Fig. [5j The scales in the figure have been adjusted to allow for the 
naive 1/M^ behaviour of the cross section, so that any departure from the 'symmetry' of the 
plot is due to the PDFs. 

At first sight, it might be surprising that there is not more difference between the predictions 
from the different sets of PDFs for such low values of X2- However, 

(a) the predictions made by using all reasonable sets of PDFs which are consistent with each 
other within the error bands; 

(b) the partons are being probed at a moderately high scale: ~ 70 and 280 GeV^ where 
the difference between different sets of PDFs is not large; 

(c) recall that, at the moment, there is no prediction for V ^ 4, but only those based on 'ad 
hoc' PDF extrapolations to very low x; however, the forthcoming LHCb data can provide 
an important direct measurement of the PDFs in this low-x region; a domain which is 
not accessible by previous data. 
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5 Conclusion 



The observation of low-mass Drell-Yan ^'^fi~ pairs at high rapidity at the LHC can probe the 
PDFs in a low x domain unreachable in previous experiments; that is, x ^ 10~^. However, the 
constraints on the PDFs found using a conventional analysis of the data, show great sensitivity 
to the choice of factorization (and renormalization) scale. Here we have introduced a technique 
to determine an optimal scale which greatly reduces this sensitivity. The optimal scale is found 
to be hf = 1.4M, where M is the mass of the produced ^~ pair. Simply for illustration, 
we show 'predictions' of the Drell-Yan cross section, as a function of rapidity, Y , for M = 6 
and M = 12 GeV for a variety of recent PDF sets, but we stress that the predictions for 
y ^ 4 depend on 'ad hoc' extrapolations. Here, LHC data can, for the first time, make a direct 
measurement of PDFs in the x 5^ 10~^ domain. 
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